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F  RE  FACE 

This  Memorandum  discusses  the  number  of  launches  required  to 
establish  and  maintain  a  spacecraft  system.  The  mean  and  dispersion 
in  number  of  launches  are  derived  for  a  system  consisting  of  an  arbi¬ 
trary  number  of  exponentially  decaying,  spacecraft  and  with  a  constant 
probability  for  successful  replenishment  launches.  The  report  shows 
the  dependence  of  these  results  on  the  spacecraft  mean  lifetime,  the 
launch  success  probability,  and  the  number  of  spacecraft  per  launch. 

The  results  arc-  useful  in  estimating  the  probable  costs  of  various 
spacecraft  systems  and  should  be  of  assistance  in  comparing  such  systems 
The  Memorandum  was  prepared  as  part  of  the  Advanced  Research 
Projects  Agency's  VELA  Analysis  study.  An  abridged  version  of  this 
Memorandum  was  published  in  the  Proceedings  of  the  IEEE:  Special 
Issue  on  Nuclear  Test  Detection,  December  1965. 
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SUMMARY 


The  mean  rc^  and  dispersion  in  nuaber  of  launches  required 
to  maintain  with  a  constant  launch  capability  a  system  of  Nq 
exponentially  decaying  spacecraft  for  a  time  t  large  compared  with 
the  mean  system  failure  time  are 


1  (  t  °1  1  m0\ 

*  PW  2m  2  2  V 


^0  P\ 

'  n1  '  V 


where  p  is  the  probability  for  a  successful  launch  and 


V'l  t 

k-N 

o 


j  -  0,1 


with  X.  the  mean  lifetime  of  the  spacecraft,  and  with  nQ  and  n^  the 
number  of  spacecraft  on  orbit  immediately  following  initial  estab¬ 
lishment  and  first  reestablishment,  respectively. 
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I.  INTRODUCTION 

In  the  course  of  some  recent  studies  of  costs  of  various  space¬ 
craft  (S/C)  syatems  for  detection  of  nuclear  explosions,  a  number 
of  useful  results  applicable  to  the  statistical  analysis  of  S/C 
replenislunent  have  been  derived.  In  general  it  is  desired  to  determine 
the  number  of  launches  required  to  maintain  a  system  of,  say,  N  S/C 
for  a  time  T  with  a  launch  capability  of  N  S/C  per  launch.  The  model 

L 

assumes  system  failure  when  the  number  of  operating  S/C  falls  below 

Nq  and  assumes  that  replenishment  launches  can  be  effected  immediately 

after  failure.  It  is  usually  assumed  that  the  probability  for  failure 

of  a  S/C  known  to  be  in  successful  operation  is  independent  of  the 

S/C  age;  hence,  the  asic  parameters  of  the  analysis  are  X,  the  S/C 

expected  lifetime  (or  MTBF) ,  and  p,  the  probability  of  a  successful 
* 

launch.  The  most  general  solution  is  the  probability  that  n 
launches  are  required  to  maintain  the  system  for  the  time  T,  although 
cost  estimates  are  usually  based  on  the  mean  and  variance  of  the 
distribution.  In  order  to  consider  possible  extensions  in  the  model, 
the  treatment  below  is  occasionally  somewhat  more  general  than  that 
required  for  the  simple  model. 

★ 

In  typical  detection  systems,  NQ  may  lie  in  the  range  1  to  10, 
Nl~1  to  5,  T~5  yr,  p  ~  0.5  to  0.9,  and  X  ~  1/2  to  3  yr. 


II.  MAINTENANCE  COMPOUND  DISTRIBUTION 


It  is  convenient  to  write  t’  .  probability  as  a  compound 
probability  distribution 


Pn<C>  '  I  £j<°  h„,j  (1) 

.1-0 

where  f^(t)  Is  t^ie  Probability  that  j  system  failures  occur  during 
the  time  t,  and  h^  is  the  probability  that  n  launches  are  required 
when  j  failures  occur.  This  approach  is  useful  since  the  events 
determined  by  the  probability  distribution  f  are  independent  of  those 
determined  by  the  probability  distribution  h. 
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III.  DISTRIBUTION  OF  NUMBER  OF  SYSTEM  FAILURES  (£.) 


The  distribution  f^(t)  for  the  number  of  system  failures  in  time 

t  will  be  derived  in  terms  of  the  probability  g(t)dt  that  the  system 

fails  in  a  small  interval  dt  at  time  t.  That  is.  g(t)  is  assumed 
* 

known,  and  f.(t)  is  to  be  determined. 

J 

It  is  convenient  to  introduce  the  probability  b^(t)dt  giving 

the  probability  that  the  j11*1  failure  occurs  in  the  interval  dt 

at  time  t.  Then,  since  the  change  df^(t)  in  f^  at  time  t  during  the 

interval  dt  is  the  difference  between  the  probability  for  entering  j 

(the  probability  that  the  j1"*1  system  failure  occurs)  during  the 

interval  and  the  prouability  for  leaving  j  (the  robability  that  the 
s  t 

(j  +  1)  failure  occurs)  during  the  interval,  it  may  be  expressed  as 


df ^(t)  =■  bj(t)dt  -  b  (t)dt 


df  ,(t) 


bJ^  "  bj+l^ 


With  the  notation  that  the  Laplace  transform  (L.T.)  of  a  function  of 

x 

time  such  as  g(t)  is  represented  by  g  (cn) ,  that  is,  if 


■w  ■  J 


dt  e  a>t  g(t) 


the  L.T.  of  Ea.  (2)  is 


Qfj*(a>)  -  b*(a>)  -  b1+1*(cn) 


In  Station  VII,  the  probability  g(t)  for  a  simple  system  is 
derived  from  the  failure  probabilities  for  the  individual  S/C. 
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It  is  useful  to  observe  that  the  L.T.  of  a  probability  density  function 
(pdf)  is  <  e  nt>,  the  expectation  value  of  e  over  the  distribution. 

The  system  is  assumed  in  operation  at  time  t  =  0  with  the  pdf 
for  failure  given  by  g^Ct).  The  interval  tc  the  first  failure  is 
represented  by  the  random  variable  Tq.  At  time  t  ■  Tq  the  system 
suffers  its  first  failure  and  is  immediately  reestablished,  after 
which  the  failure  pdf  is  g^(t).  The  process  continues,  with  intervals 
between  failures  represented  by  the  independent  random  variables 
T0 »r i ,  •  •  • »  with  failures  at  times  t^t^,...  given  by 

C  j  ^  Tk 

k»0 


and  with  giving  the  failure  pdf  appropriate  to  the  interval 

following  the  k^  reestablishment.  It  is  convenient  to  include  the 
j  “  0  term  with  the  assumption  t^  •  0.  Then,  since 


it  follows  that 


and,  for  j  £  1 


V“> 


b^  (cu) 


-  1 


The  functions  g(t)  and  b^(t)  give  probabilities  p“r  unit  time  and 
are  density  functions,  whereas  f^(t)  and  Pn(t)  give  probabilities  in  i 
and  n  and  are  not  density  functions. 
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where  the  last  step  follows  from  the  Independence  of  the  intervals 


r,  .  Thus 
k 


b  j  (a:) 


r 

1 

j  =  o 

1  j_1 

TT  gk*(03) 

j  *  1 

k=0 

and,  from  Eq.  (4) 


fjV) 


1 


★  1 

8,(03)-,  J*  . 

L  ^  J  ^  8k  (®) 

aig  .  (oj)  k«0 


(5) 


Equation  (5)  gives  the  L.T.  of  the  failure  distribution  in  terms 
of  the  failure  pdf's  g^(t). 
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IV.  DISTRIBUTION  OF  NUMBER  OF  LAUNCHES  FOR  FIXED  NUMBER  OF  FAILURES  (h  4) 

_ _ _ .  I. . n»J_ 


uk 


If  the  system  fails  exactly  j  times,  and  if  the  random  variable 

represents  the  number  of  launches  required  to  reestablish  the 
th 


system  after  the  k  failure,  then  h  .  is  the  probability  that  S  “  n, 

n » J  j 

where  is  given  by  the  sum 


r  I 


(6) 


k“l 

Consider  first  the  sum  of  two  non-negative  integral  random  variables 
and  If  an£^  represent  respectively  the  probabilities 

that  «  i  and  *  j,  then  the  probability  that  the  sum  +  u2' 
equals  n  is  given  by 


h  -  V  41)  «<2> 

n  Lt  k  n-k 
k"0 

Let  A^  \s)  be  the  probability  generating  function  (pgf)  for  a^k\ 


(7) 


defined  by 


A^k\s)  *  ^  a^Ck)^ 


(8) 


i»0 

Then,  from  Eq.  (7),  the  generating  function  H(s)  for  h^  is  given  by 

H(s)  -  A^'(s)  A^2)(s) 


The  generalization  of  this  result  to  higher-order  sums  is  obvious: 

(k) 

if  a^  is  the  probability  that  i  launches  are  required  to  reestablish 


See  Ref.  1,  p.  251. 


f  •-» 


-  -r  - - 
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the  system  after  the  kth  failure,  the  generating  function  for  the 
distribution  in  total  number  of  launches  required  to  maintain  the 
system  through  j  failures  is  given  by  the  product 

j 

H  (s)  *  IT  A(k)(s) 

J  k-0 

For  completeness ,  the  terra  for  j  =  0  is  included;  after  the  "zeroth" 
failure,  zero  launches  are  required  to  reestablish,  that  is, 

=  0  for  all  i  except  i  *  0,  whereas  a^^*  A^(s)  =>  Hq(s)  -  1 


(9) 
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V. 


GENERAL  SOLUTION  FOR  DISTRIBUTION  OF 
LAUNCHES 


TOTAL  NUMBER  OF  REPLENIS HME NT 


In  Section  III,  the  L.T.  of  f^(t),  the  number  of  system  failures 
in  time  t,  was  determined  in  terms  of  the  L.T.  of  the  system  failure 
pdf  s  g^(L);  and  in  Section  IV,  the  pgf  for  the  probabilities  h 

n»  j’ 

the  distribution  in  number  of  launches  for  a  fixed  number  of  system 

failures,  was  determined  in  terms  of  the  pgf  for  the  probabilities 

(k)  .  .. 

ai  ,  the  distribution  in  number  of  launches  required  to  reestablish 
after  a  single  system  failure.  From  Eq.  (1),  the  L.T.  of  the  pgf  for 

Pn^’  the  total  number  of  replenishment  launches  required  in  time  t, 
is  given  by 

CO 

★  P  ★ 

p  (tD,s)  -  )  fj  (o>)  Hj(s) 

J=0 

Hence,  from  Eqs.  (5)  and  (9) 


P*telS)  .  £  [ililiii]  m  7T  AU>(S) 

j=o  i  (co)  k"°  i-0 


(10) 


Since  the  functions  g^  (to)  and  A^k\s)  are  given  in  Eqs.  (3)  and  (8), 

in  terms  of  known  functions  g^(t)  and  respectively,  Eq.  (10) 

respresents  the  general  solution  for  the  distribution  P  . 

n 


t 


/tar” - 
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VI.  TOTAL  REPLENISHMENT  DISTRIBITCI.ON  FOR  CONSTANT  LAUNCH  SUCCESS 
PROBABILITY 


The  model  described  In  the  Introduction  assumes  a  constant 
probability  for  a  successful  launch.  Hence,  the  replenishment  launch 

(k) 

number  probabilities  do  not  depend  on  k,  the  system  failure 

number,  and  may  be  written  as 


(k) 


k  £  1 


where  a^i  the  probability  that  i  launches  are  required  to  replenish 
the  system  after  failure  (that  is,  the  probability  that  the  i1"*1 
launch  is  the  first  successful  launch  following  a  system  failure)  is 
given  by 


a  =  < 
i 


pq 


i-1 


i  =  0 

i  2  1 


(11) 


and  q  ■  1  -  p.  The  pgf  for  the  probabilities  a^  is 


00 


A(s)  -  p  £  q1_V  -  (12) 

i**l 
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Expansion  of  Hj(s)  in  a  power  series  in  s  gives 
Hj(s)  »  (ps)1  £  — +-  k  '  1)J 


k»0 


(j-1).'  k.' 


(qs)' 


, — 1 


n**  j 


InzllL 


(j-1)'  <n-j) - 


nn*J  nJ 
q  p 


(14) 


Hence,  the  probability  that  n  replenishment  launches  are  required 

if  j  failures  (j  *  1)  occur  is,  from  Eq.  (U) 


r 


n,  j 


n  <  j 


Cn-1).' 


(j-1)'  (n-j).'  q  J  pJ  n  *  j 


If  j  =  0  (no  failures),  no  replenishment  launches  are  required; 

therefore,  h^  ^  =  1  and  hn  ^  «  0  for  n  £  1.  Substitution  of  this 

expression  for  h  into  Eq.  (1)  gives 
n ,  j 


(15) 


f0(O 


Pn(t)  «<!  n 


\  (n-1) n-j  j  ,  . 

L  (j-i>:  (ii- j).'  q  p  fj(0 


n  =  0 


n  S  1 


j-1 


(16) 
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VII.  AGELESS  SPACECRAFT 


Consider  a  system  consisting  of  a  single  S/C  placed  in  operation 
at  time  t  *  0,  and  let  ^(t)  represent  the  probability  that  it  remains 
in  successful  operation  (is  "alive")  at  age  t.  Assume  that  the  S;C 
is  observed  to  be  alive  at  age  tQ,  i.e.,  ^ ( t )  =  1  for  0  £  t  £  tQ. 

For  times  beyond  the  last  observation  time  t^,  the  decrease  di)  in  i(t) 
in  an  interval  dt  ^the  probability  that  the  S/C  fails  in  dt)  is  the 
product  of  the  probability  that  the  S/C  is  alive  at  age  t  and  the 
conditional  probability  ft)(t)dt  that  the  S/C,  if  alive  at  age  t,  fails 
in  the  interval  dt 

-di  *  i'(t)  cp(t)dt 

With  the  given  boundary  conditions,  ^(t)  becomes 


*(t) 


fl  0  *  t  S 

|  -r  dt  <KO 


t 


0 


For  convenience,  and  because  of  reasonable  agreement  with 
experience,  it  is  usually  assumed  that  the  failure  probability  of  an 
operating  S/C  is  independent  of  its  age;  hence 


<P(0 


1 

X 


where  X  is  a  constant,  and 


^(t)  « 


07) 


from  which  it  is  easily  determined  that  X  is  the  MIBF  measured  from  t^ , 
the  (most  recent)  time  of  known  successful  operation. 
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In  general,  the  probability  'i'j(t)  that  t  system  of  several  S/C 
remains  in  successful  operation  for  some  time  t  following  its 
reestablishment  is  one  minus  the  probability  that  it  has  failed 


during  the  interval  t 


i.  e. 


With  the  L.T.  of  this  relation 


(18) 


Yj*((o) 


may,  from  Eq.  (5),  be  expresseu  as 


j  *■  0 

J  *  1 


(19) 


(20) 


Thus,  fj(t),  the  probability  that  there  are  exactly  j  failures  in 
time  t,  is  the  convolution  of  the  probabilities  for  j  failures  g^ , 


g^,....g  ana  the  probability  Y 
jth  failure. 


for  continued  operation  after  the 


The  model  assumes  a  constant  launch  capability  and  a  constant 
orbit  requirement  Nq.  Consider,  for  example,  a  system  with  Nq  ■  5 
and  N  =  3.  Two  successful  launches  ars  required  to  establish  the 

L 

system;  vence,  at  t  =  3,  there  will  be  six  operating  S/C.  At  the 
first  system  failure  (i.e.,  at  the  second  S/C  failure),  one  successful 


launch  is  needed  to  replenish,  immediately  after  which  there  will  be 
seven  operating  S/C.  Successive  failures  form  a  recurrent  pattern 
between  four  and  seven  S/C  on  orbit. 


In  general,  there  will  be  n^  S/C  in  operation  after  initial 
establishment  and  n^  S/C  immediately  after  the  first  and  all  subsequent 
replenishments.  If,  in  addition,  the  S/C  are  ageless,  then  all  system 
failure  pdf's  after  the  first  failure  are  identical,  i.e. 


8j(t) 


*0(t) 


j  -  0 


j  *  1 


(21) 


With  this  result, 


Eq.  (20)  becomes 


8, 


0 


r 

(8l  ) 


_  1 


j  -  o 

j  *  1 


(22) 


Following  initial  establishment,  the  number  of  S/C  in  successful 

operation  decreases  from  nQ  through  successive  S/C  failures.  However, 

the  system  is  considered  to  fail  only  when  the  number  of  S/C  falls 

below  Nq.  Hence,  the  probability  YQ  for  successful  system  operation  is 

the  sum  of  the  probabilities  that  k  S/C,  where  N  Ski  n„,  are  in 

o  0 

operation;  i.e. 


I 

k*»N 


n„ 


Vk 


(23) 


and  Y^  is  given  by  the  same  expression  with  n^  replaced  by  n^.  In 

the  following,  the  subscripts  0  and  1  will  often  be  dropped  from  Y, 

g,n  and  associated  parameters  when  the  results  apply  to  either 

subscript.  The  L.T.  of  Eq.  (23)  depends  on  the  integral 

00  .  L  k  .  t  n-k 

f  e-cut:(e  X)fle  X)  dt 

J0  ,  /  \  / 

which,  with  the  substitution  x  ■  e  ,  may  be  transformed  into  the 


integral 


x  J  xa*-+k'1n-x)r,-k  dx 


This  integral  defines  the  beta  function^  and  therefore  may  be 


expressed  as 


r<m\+k)  T(n-k+l) 
r(a&+n+l) 


where  T  is  the  gamma  function.  Thus,  the  L.T.  of  Y  is 


Y  (o>)  -  X 


nl _  V  r 

r(ojx-Hi+i)  L  k.‘ 


Except  in  the  case  nQ  -  n^  *  N  ,  inversion  of  the  expression  Eq.  (22) 
,  .  *  * 

for  t  with  the  above  value  of  Y  is  not  feasible,  hence,  the 

replenishment  distribution  P  cannot  be  determined  in  detail,  and  it 

n 

is  necessary  to  evaluate  its  mean  and  variance  by  other  methods. 
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VIII.  MEAN  AND  VARIANCE  OF  REPLENISHMENT  DISTRIBUTION 


From  Eq.  (1),  the  mean  of  the  replenishment  distribution  Pn 


is  given  by 


”R=  I  "  Pn-  I"  l  f3  hn.j 

n“0  n“0  j=0 


■  I  ■■  <s>> 


Where  <^Sj^>,  given  by 


^S3>“  l  "  Vj 


is  the  expected  total  number  of  launches  required  when  j  failures 


occur. 


Thus,  <^S^  ^  is  the  expected  value  of  the  sum  S^  given  by 


Eq.  (6),  i.e. 


/  j  \  -i 

<sj)  ■(!  \)  ■  I<"r>* 

'k*<l  /  k-1 


where  m&  is  the  mean  number  of  launches  required  to  reestablish  a  system 
after  a  single  failure.  The  last  step  in  Eq.  (26)  follows  from  the 
assumption  that  the  distributions  in  number  of  reestablishment  launches 
are  independent  of  k.  Thus,  from  Eqs.  (25)  and  (26) 


"t  *  I  f]  h  ■  "a  ’f 


(27) 


1.6 


where 

in  -  «  }  j  1, 

£  L  J 

3*o 

is  the  mean  number  of  system  failures.  A  similar  calculation  for 
the  variance  in  total  number  oi  replenishment  launches,  defined  by 


★ 

gives 


:r> 


n:*0 


m.a 

f  a 


+  in 


2J 

a  f 


(28) 


where  is  the  variance  in  the  number  of  launches  required  to 

2 

reestablish  after  a  single  failure,  and  is  the  variance  in  number 
of  system  failures.  The  distribution  in  number  of  launches  to 
reestablish  after  a  single  failure  is  given  by  a^  in  Eq.  (11);  the 
mean  and  variance  of  this  distribution  are 


m  = 
a 


1 

P 


(29) 


2 


c 

a 


S 

0 

P 


(30) 


where  p  is  the  probability  of  a  successful  launch  and  q  =  1-p. 
The  expressions  Eqs.  (27) -(30)  may  be  combined  to  give 

mf 

and 

2  CTc  + 

O'  =  — - 

R  2 

P 


(31) 


* 


See  Ref. 


1,  p.  276 


(32) 
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From  Eq.  (22),  the  L.T.  of  the  mean  m^  is 


■  <»*  ■  i  *; 

J-0 


*  *  r*  * 

\  80  l  J<«!  > 

j-l 


j-l 


yu  *  a  *  __d_  V  /  O 

1  0  dR  *  l  (81  } 

81  j“l 


Y  *  g  *  — sL- 

1  80  * 


dg 


1 


[d-S1*>'1-  1] 


,  *  *)‘2  g0 

'  <1  g0  U-*!  "  17 

0) 

2 

The  variance  a  may  be  obtained  from  the  expected  value  of  J(J+1); 
the  L.T.  of  this  quantity  is 


(33) 


<J(j+D>*  -  V*  gQ*  £  J(j+l)(g1*) 

J-l 


j-l 


«  y  *  g  *  -i 

1  80 


dg  *2  4, 

81  j- 


I  <0 


J+l 


V  g0*  [S^CI-S,*)-1-  h*] 

dg.  u  J 


„  *  ★  *  -3 

2Y1  80  (1'81  )  "  3.  *  2 

CD  (^  ) 


(34) 


18 


Inversion  of  Eqs.  (33)  and  (34)  with  the  relation  Eq.  (24)  for 

★ 

Y  is  not  feasible;  however,  approximate  expressions  for  the  inverse 
function  may  be  obtained  in  the  most  important  limit,  that  of  large 
times,  with  term  by  term  inversion  of  the  expansion  of  the  L.T.  in  a 
power  series  in  uu. 

Repeated  application  of  the  recursion  relation  for  the  gamma 
function,  i.e. 

T(l+x)  =  x  r(x) 

gives 


r(kdtcX)  =  (k-l+trX)  (k-24tuX)  (1+tnX)  ^(l-taX) 

k-1 

-  (k-1)!  +  aiX(k-l) !  j  +  0(te2)]  r(l-ttox) 

j=l 


Also,  since 


r(i4u>x)  =  r(i)  +  mx  r  '(i)  +  o^2) 


it  follows  that 


(k-Ty  =  1  +  to(ak+c)  +  0(cu2) 


where  c  is  a  constant  and 


k-1 


>:  i 

J-l 


j 

* .  2 


Then,  from  Eq.  (24),  'f  (<jd)  is  given  to  order  £jd  as 


Am) 


l-rrn(0f  ,+c)  ^ 

k»N 


l-+tD(Ofk+c) 


x  I 

k=N 


l-cD(an+ro,k) 


k 


m  -  Yl- 


(35) 
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where 


m  =  X 


u 


Y  “  X  ) 


n  n  n 

r  an+i"ak  .  2  r  i  n 


lUl) 

k«N  j*k 
o  J 


The  significance  of  the  parameters  m  and  y  (i.e.,  of  the  four  parameters 

m^.m^jY^,  and  Y2)  may  be  seen  from  a  comparison  of  the  two  previous 
★ 

expressions  for  g 


★  ★ 
g  =  1  ■  a  I 


To  order  03  ,  these  relations  are  respectively 


*  2 
g  *  1  -  can  +  ya> 


Thus, 


1  -  «.<t>  +  a-  <t2> 


<t> 


is  the  mean  system  life,  and 


<t2> 


The  dispersion  a  in  system  lifetime  is 
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<t2>-<t>2-  2V- 


m 


n  n 


k-N  j-k 
o  J 


k»N  j=N 
o  o 


k-N 


*  - 1 
With  these  results,  m,  may  be  expanded  to  order  uu  as 


8 


m. 


0 


1-cun 


0 


2  * 

cd  ^  2  , 

1  CD  (.»  -V,CD) 
▲  L 


-r —  +  -  (— j  -  — ) 

2  CD  V  2  m,/ 

cd  m1  Tn^  1 


(39) 


Insertion  of  this  expression  gives  the  mean  number  of  failures  in 
time  t  to  order  constant,  that  is,  neglecting  terms  of  order  m^/t 
and  ra^/ t  as 


t  CT1  1  mo 

m  -  + - r  +  - - 

1  “l  2T.J2  2  ml 


Similarly,  to  order  cd 

<ju+d> 


26, 


0 


_ ~<gnQ) 

CD3('i'1*)2  CD3(m1-V1CD)2 

2V 


3  2 

UTn^ 


[1+  a(_J  .  mQ)J 


2  2  _1_  V 

3  2  +  2  (1  +  2  ~  m 

(D  m  1 


(40) 


i 
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Inversion  gives,  to  order  t/m 


<JOD>  -  -^  +  ^U+- 


2  m. 

tn^  1 


To  order  t/m,  the  variance  in  number  of  failures  is 


t  r  m0.“l 

-  ~  1  +  2(~-  -  —) I 

ml  L  m  2  ml  1 


These  results,  together  with  the  results  shown  in  Eqs.  (29)  and  (30), 
may  be  used  in  Eqs.  (27)  and  (28)  to  give  the  mean  and  variance  in 
total  number  of  launches  required  in  the  (large)  time  t  as 


i  (x  +  JT.  +  A  . 

»  X  a,2  2  "i 


The  expression  Eq.  (42)  is  accurate  to  the  constant  term;  the 
expression  Eq.  (43)  is  accurate  to  order  t/m. 

For  the  important  special  case  ■  1,  the  distribution  f ^ 
may  be  obtained  exactly.  With  one  S/C  per  launch,  there  will  be 
no  surplus  S/C  in  either  the  initial  or  reestablished  phases,  i.e., 
n^  =  n^  =  Nq.  From  Eqs.  (36)  and  (37)  the  mean  and  variance  in  the 
system  lifetime  following  establishment  (and  each  subsequent 
reestablishment)  are 


and 


From  Eq.  (24) 


and,  from  Eq.  (22) 


2  mK_ 

N 


Ao>)  x 


m 


m 


cd^+N  airrfl 
o 


*  *  J 

f  -  *  (l-a.'J  )  -  -  , 

1  (ran+l)J 


m 


Inversion  gives  the  Poisson  distribution 


t 

m 


j 


e  ,JL 

11  V 


2  *>  2 

Note  that  when  the  relations  =  mg"  -  m 


^  ,  which  foil  tw  from  the 
assumption  -  1,  are  used  in  the  approximations  Eqs.  (4C)  and  (41) 
for  the  mean  and  variance  in  £,  the  results  are 


m. 


m 


in  agreement  with  the  mean  and  variance  of  the  exact  (Poissoi) 


distribution 
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IX.  MEAN  AND  VARIANCE  OF  ESTABLISHMENT  AND  TOTAL  DISTRIBUTIONS 


The  probability  distribution  in  number  of  launches  required  to 

ertablish  a  system  is  given  essentially  by  h  the  probability 

n»  J 

that  n  replenishmen  launches  are  required  if  j  system  failures  occur. 

The  establishment  distribution  is  obtained  from  t^e  expression  for 

h^  given  in  Eq.  (13)  if  'he  random  variable  n  represents  the  number 

of  establishment  launches  and  the  integer  j  represents  the  minimum 

number  of  successful  launches  required  to  establish  the  system.  The 

2 

mean  in  and  variance  a  in  number  of  establishment  launches  are, 

E  E 

according  to  Eq.  (15),  given  bv 


m 


E  P 


(44) 


and 


2  _ o  o 

aE  "  J.  ~  ~ 


(45) 


where  j  has  been  replaced  by  Lq. 

If  N  represents  the  total  (establishment  plus  replenishment) 

number  of  launches  required  in  time  t,  the  maan  m^  and  dispersion 
2 

axl  in  N  are  giver,  by 
N 


“n  ■  mE  +  "r 


(46) 


and 


2  2  ,  2 
aN  '  aE  +  JR 


(47) 


where  and  are  given  by  Eqs.  (42)  and  (43) . 


DOCUMENT  CONTROL  DATA 


I  ORK3NAT1NC  activity 

the  RAND  CORPORATION 


2o  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2b  GROUP 


3  REPORT  TITLE 


STATISTICAL  ANALYSIS  OF  SPACECRAFT  REPLENISHMENT 


4  AUTHOR! S)  (Uxl  rwm* ,  firtt  nonw.inltlol) 

Modesitt,  G.  E. 


3  REPORT  DATE 

March  1966 


7  CONTRACT  or  GRANT  NO. 

SD-79 


to.  TOTAL  NO.  OF  PAGES 

37 


8.  ORIOINATOR  S  REPORT  NO. 

RM-4739-ARPA 


6b  NC.  OF  REFS 

2 


l  W  )  1 B"  1  lhi7i  R  r  U  f  i  <  M  4- 


3.  ABSTRACT 

A  discussion  of  the  number  of  launches 
required  to  establish  and  maintain  a 
spacecraft  system,  undertaken  to  estimate 
the  probable  costs  of  various  spacecraft 
systems  and  to  be  of  assistance  in 
comparing  them. 


9b  SPONSORING  AGENCY 

Advanced  Research  Projects  Agency 


II  KEY  WORDS 

Spacecraft 
rlaintenance 
Launching  systems 
Interplanetary  flight 
Space  technology 
Statistical  methods  and 
Procedures 


